Ground-State and Domain- Wall Energies in the Spin-Glass Region 
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The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond 
Ising model on square lattices with independent, identically distributed bonds of probability p of 
Jij = — 1 and (1 — p) of Jij = -1-1 are studied. We are able to consider large samples of up to 320^ 
spins by using sophisticated matching algorithms. We study L x L systems, but we also consider 
LxM samples, for different aspect ratios R = L/M . We find that the scaling behavior of the ground- 
state energy and its sample-to-sample fluctuations inside the spin-glass region {pc < p < 1 — Pc) 
are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusp-like 
singularity at pc. Inside the spin- glass region the average domain- wall energy converges to a finite 
nonzero value as the sample size becomes infinite, holding R fixed. Here, large flnite-size effects 
are visible, which can be explained for all p by a single exponent u ~ 2/3, provided higher-order 
corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for ii — > 0: 
the distribution of the domain-wall energies converges to a Gaussian for i? 0, although the domain 
walls of neighboring subsystems of size L x L are not independent. 
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I. INTRODUCTION 

The spin glass (SG) phase^'^ is not stable at finite tem- 
perature in two dimensions (2D). When the energies are 
not quantized, the behavior at T = can be understood 
in terms of a scaling theoryt^i^ which is usually called the 
droplet model. This theory describes the scaling of the 
average domain-wall (DW) energy Edw with length scale 
in terms of the stiffness exponent 6, i.e. 



\Edv 



(1) 



This exponent has a value of about —0.28 for 2D, al- 
most independent of the detailed nature of the bond 
distribution^ and describes the scaling of different kinds 
of excitations like domain walls and droplets, at least 
if large enough system sizes are studied.^ The possibility 
that quantization of the energies might lead to the special 
behavior 9 = was first pointed out by Bray and Moore.— 
It has become clear over the last few years that there ac- 
tually is a fundamental difference in the behavior of the 
2D Ising spin glass at zero temperature, between those 
cases where the energies are quantized and those where it 
is notj^i^i^ Nevertheless, recent result o^°i^^ indicate that 
in the low-temperature critical scaling regime, the be- 
havior for quantized and non-quantized models might be 
very similar. 

The Hamiltonian of the Edwards- Anderson model for 
Ising spins is 



H 



(2) 



where each spin Si is a dynamical variable which has two 
allowed states, -fl and —1. The (ij) indicates a sum 



over nearest neighbors on a simple square lattice of size 
L X M. The standard model for quantized energies is the 
± J model, where we choose each bond Jij to be an in- 
dependent identically distributed (iid) quenched random 
variable, with the probability distribution 

P(J,,) = pS{J,j + 1) + (1 - p)SiJ,j - 1) . (3) 

Thus we actually set J = 1, as usual. The concentration 
of antiferromagnetic bonds is p, and (l—p) is the concen- 
tration of ferromagnetic bonds. With the P{Jij) of Eqn. 
([3]), the EA Hamiltonian is equivalent to the Z2 gauge 
glass model/i^ Wang, Harrington and Preskill^ have ar- 
gued that the anomalous behavior of the ± J model is 
caused by topological long-range order, as a consequence 
of the gauge symmetry. 

Along the T = axis this model exhibits a phase 
transition from a ferromagneticly ordered phase'^^ for 
small concentrations p < Pc oi the antiferromagnetic 
bonds to a SG critical line at large values p > Pc 
(and p < 1 — Pc due to the bipartite symmetry of the 
square lattice). Recently, this phase transition was char- 
acterized by high-precision ground-state calculations ^iii^ 
which have in particular yielded Pc = 0.103(1). The fer- 
romagnetic phase persists at finite temperatures T > 
for p < pc{T), while the spin-glass correlations become 
long-range only at zero temperaturei^i^ii^ Interestingly, 
Pc{T) > pc(0)r^i^'^ i.e. the ferromagnetic phase is reen- 
trant. 

The first hint of the remarkable behavior of the ±J 
model in the SG regime was observed by Wang and 
Swendsen.-^ They found that, although for periodic 
boundary conditions there is an energy gap of 4 J between 
the ground states (GS) and the first excited states, the 
specific heat when T/J <C 1 appeared to be proportional 
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to exp(-J/2r). 

This result was questioned by Saul and Kardar. ^^i^^ 
The behavior of the specific heat was finally demon- 
strated in a convincing fashion by Lukic et alP'^ Saul 
and Kardar also found that the scaling of the DW en- 
tropy with size did not appear to agree with the predic- 
tion of the droplet model. This issue has recently been 
clarified, and the DW entropy scaling anomaly has been 
associated with zero-energy domain walls. 

In this work we explore the behavior of the ± J model 
inside the full SG phase Pc < P < 0.5 (the behavior 
0.5 < p < 1 — is equivalent due to the symmetry of 
the model). In particular we study the finite-size scaling 
behavior of the GS energy, of the fiuctuations of the GS 
energy, and of the domain-wall energy. We also employ 
aspect-ratio (AR) scaling, ^■^'^^ i.e. we study rectangular 
lattices of width L and height Af , the aspect ratio being 
R = L/M. Carter, Bray and Moorc^^ have extended AR 
scaling to the Ising SG, and demonstrated that study- 
ing the scaling as a function of R is an effective method 
for calculating the exponent 9. Here, we look at the limit 
i? ^ and show that AR indeed works for the ± J model 
as well, in contrast to previous attemptsi^ The entire 
probability distribution P{Edw) has a simple Gaussian 
form for small R, although, as we will show, the underly- 
ing assumption of independent contributions to the DW 
energy is not strictly valid. 



II. METHODS 

We define domain walls for the SG as it was done in 
the seminal work of McMillan. We look at differences 
in the GS energy between two samples with the same 
set of bonds, and the same boundary conditions in one 
direction, but different boundary conditions in the other 
direction. In particular we use periodic (p) or antiperi- 
odic (ap) boundary conditions along the x-axis and free 
boundary conditions along the y-axis. For an example of 
a DW created in such a way, cf. Fig. [TT] For each set of 
bonds, the DW energy Edw is then defined to be the dif- 
ference in the GS energies of the two different boundary 
conditions: 



Ed 



En ~ En 



(4) 



We use a so-called matching algorithm to calculate the 
Qg ^27,28 sketch just the basic ideas of the 

matching algorithm. For the details, see Refs. [27ll29l . 
I30II31I . The algorithm allows us to find ground states 
for lattices which are planar graphs. This is the reason 
why we apply (anti-)periodic boundary conditions only 
in one direction, x, while the other direction, y, has free 
boundary conditions. In the left part of Fig. [1] a small 
2D system with (for simplicity) free boundary conditions 
in both directions is shown. All spins are assumed to 
be "up" , hence all antiferromagnetic bonds are not sat- 
isfied. If one draws a dashed line perpendicular to each 
broken bond, one ends up with the situation shown in 
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FIG. 1: (color online) 2D Ising spin glass with all spins up 
(left, up spins not shown). Straight lines are ferromagnetic, 
jagged lines are anti-ferromagnetic bonds. The dashed lines 
connect frustrated plaquettes (crosses). The bonds crossed 
by the dashed lines are unsatisfied. In the right part the 
GS with three spins pointing down (all others up) is shown, 
corresponding to the minimum number of unsatisfied bonds. 



the figure: all dashed lines start or end at frustrated 
plaquettes and each frustrated plaquette is connected to 
exactly one other frustrated plaquette by a dashed line. 
Each pair of plaquettes is then said to be matched. In 
general, closed loops of broken bonds unrelated to frus- 
trated plaquettes can also appear, but this is possible 
only for excited states. Now, one can consider the frus- 
trated plaquettes as the vertices and all possible pairs of 
connections as the edges of a (dual) graph. The dashed 
lines are selected from the edges connecting the vertices 
and called a perfect matching, since all plaquettes are 
matched. One can assign weights to the edges in the 
dual graph, the weights are equal to the sum of the abso- 
lute values of the bonds crossed by the dashed lines. The 
weight A of the matching is defined as the sum of the 
weights of the edges contained in the matching. As we 
have seen, A counts the broken bonds, hence, the energy 
of the configuration is given by 



(5) 



{iJ) 



with \Jij\ = 1 from Eqn. Note that this holds for 

any configuration of the spins, if one also includes closed 
loops in A, since a corresponding matching always exists. 

Obtaining a GS means minimizing the total weight 
of the broken bonds (see right panel of Fig. [T]). This 
automatically forbids closed loops of broken bonds, so 
one is looking for a minimum-weight perfect matching. 
This problem is solvable in polynomial time. The al- 
gorithms for minimum- weight perfect matchings^i^^ are 
among the most complicated algorithms for polynomial 
problems. Fortunately the LEDA library offers a very 
efficient implementation^ which we have applied here. 

With free boundary conditions in at least one direction, 
Edw is a multiple of 2 (since J = 1). When periodic or 
antiperiodic boundary conditions are applied in both the 
X and y directions, however, Edw becomes a multiple 
of 4 if the number of —1 bonds is even. If the number 
of —1 bonds is odd, Edw takes on values (4n -I- 2) with 
these boundary conditions. When M is odd, changing 
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the boundary conditions in the x direction from periodic 
to antiperiodic changes the number of —1 bonds from 
odd to even, or vice versa. 

The behavior of the DW energy for the ±J model 
with R > 1 was discussed previouslyjSS For the bound- 
ary conditions we are using, the average (iSdiul) goes to 
zero exponentiaUy as R becomes much greater than one. 
Here (...) denotes an average over the ensemble of ran- 
dom bond distributions for an L x M lattice. We can 
understand this result by thinking of the system as con- 
sisting of blocks of size M x M pasted together along 
the X direction. The probability of having a zero-energy 
domain wall in each block is almost independent of the 
other blocks. The same thing happens when M is even if 
the boundary conditions in the y direction are periodic 
or antiperiodic. However, if M is odd and the boundary 
conditions in the y direction are periodic or antiperiodic, 
{\Edw\) goes to 2 at large R, because then Edw = is not 
allowed. Since a critical exponent should be independent 
of boundary conditions, this is a(nother) demonstration 
that 6* = for the ± J model in 2D. 

As pointed out earlier—, these results for large R do 
not agree with the scaling law prediction of Carter, Bray 
and Moore^^ due to the special role of zero-energy do- 
main walls. On the other hand, in the limit R —^ with 
our boundary conditions, the prediction of Carter, Bray 
and Moore is 



samples used varies with size, ranging between typically 
100000 (L = 4) to typically 10000 {L = 320). 



(6) 



This limit has not been studied before for the ± J model. 
When e ^ 0, {\Edw\) should scale as R'^/^ in 2D. For 
small R we may think of the lattice as consisting oi LxL- 
sized subsystems stacked in the y direction, with Edw 
being the sum of the DW energies of the subsystems. 
Therefore, by the central limit theorem, we anticipate 
that the probability distribution of Edw should approach 
a Gaussian distribution in the limit of small R. In the 
spin-glass region of the phase diagram, the center of this 
limiting Gaussian will approach zero as L increases. We 
will show that our numerical results for the small R limit 
are indeed in good agreement with these expectations, 
even though the subsystems are not fully independent. 



III. NUMERICAL RESULTS 

A. Ground-state energy 

We begin the study of the T = behavior of 
the ±J random bond model by studying the GS 
energy for different concentrations p of the antiferro- 
magnetic bonds [p = 0.05 (ferromagnetic phase), p — 
0.11, 0.12, 0.13, 0.14, 0.15, 0.17, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45 
and p = 0.5] and different system sizes L — 4, . . . , 320 
for square samples N = L x L, i.e. aspect ratio R = I. 
All results are averages over many different realizations 
of the disorder. The minimum number of independent 




1000 



FIG. 2: (color online) The symbols show the average ground 
state energy eo{L) for the square systems {R = 1) as function 
of the system size L for selected values of p. The lines give the 
results to the fits {L > 16) (see text). The data for p = 0.05 
look similar, but is outside the frame (eo = —1.8025). The 
inset shows the same data, rescaled as (eo(L) — eo) x L so 
that the expected behavior yields straight lines when plotted 
as function of 1/L. The data for p = 0.05 have been shifted 
downwards by 0.2 for better visibility. 

In Fig. O the GS energy per spin eo(L) is shown as a 
function of the system size for selected values of p. This 
finite-size scaling behavior yields particular insight into 
the physics and seems to be related^i^ to the stiffness 
exponent 9. It has been argued^ that for the case of 
open boundary conditions in exactly one direction and 
periodic boundary conditions in the other direction the 
GS energy per spin follows to lowest orders the form 



bed 
eQ{L)^eQ + - + — + — 



(7) 



Note that the arguments used in Ref. |35| should apply 
for the ferromagnetic phase p < 0.103 as well. Since we 
expect 9 = everywhere in the SG phase (see below) 
and 9=1 for the ferromagnetic phase, we can restrict 
ourselves to the contributions eo, b/L and c/L"^. When 
fitting^ the data to this eo{L), fits with very high qual- 
ity result, as shown by lines in Fig. [2] This is confirmed 
when plotting the data and the fitting functions in the 
form {eo{L) — eo)L as a function of 1/L. This should give, 
according to Eqn. ([7]), a linear behavior with slope c and 
ordinate intersection b. We see that indeed the data fol- 
lows the functional form well, although even higher order 
corrections seem to be present. They become significant 
at very small sizes, but we cannot quantify them within 
the statistical accuracy of the data. We only observe that 
this correction term seems to have an opposite sign for 
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FIG. 3: (color online) The symbols show the width a of the 
distributions of the GS energy eo for square systems {R = 1) 
as function of system size L for selected values of p > Pc- The 
data is rescaled with the average GS energy eo(L) and with 
the leading behavior L. Additionally, the dashed line displays 
the result for p = 0.05. The full lines show fits to constants 
at larger system sizes. 



small {pc < p < 0.25) compared to larger values of p. 
Furthermore, the behavior of the c/L^ term is quite sim- 
ilar inside the SG phase for different values of p, while 
the behavior in the ferromagnetic phase [p = 0.05) is 
very different. The value of the ordinate intercept h is 
monotonic in p, as expectedi^ 

Next, we consider the width a of the distribution of 
GS energies per spin. For short-range d-dimensional spin 
glasses, it has been proven^l that a ~ L^"^/"^, i.e. a ^ 
in this case. To remove finite-size effects caused by 
our special choice of boundary conditions, we normalize 
a by the finite-size GS energy. Thus, we plot La/\eo{L)\ 
as a function of system size and expect a horizontal line 
at Ca{p) = limL^aa La / eo{L) . As seen in Fig. [3l this is 
indeed the case. Only corrections for very small system- 
sizes become visible. We have also included the data 
for p = 0.05, which behave in the same way. Hence, 
here the ferromagnetic phase and the SG phase cannot 
be distinguished just from the asymptotic behavior, in 
contrast to the behavior of the mean alone. 

On the other hand, as visible in Fig. [31 the value of 
Ca-{p) increases when approaching the phase transition 
p Pc- For a more detailed analysis, we have plotted 
Ccr{p) as function of p inside the SG phase in Fig. 21 also 
for values of p closer to Pc than those shown in Fig. [31 One 
can see a strong increase when approaching the critical 
concentration. Fitting a power law (with pc = 0.103) 



c„(p)^c, + Kip-p,y 



(8) 



in the range p e [0.1 . 
K = -0.38(4) and n 



.0.17] yields values = 0.422(3), 
= 0.64(5), i.e. a cusp at the phase 




FIG. 4: The symbols show the limiting value Co-(p) = 
limL^cx) La/eo{L) for the width a of the distributions of the 
GS energy (as shown in Fig. [Si as a function of p. The solid 
line represents a fit according to Eqn. see text. 



transition. Note that the value of the scaling exponent is 
very close to the exponent uj = 2/3, which describes the 
finite-size scaling of the DW energy -Ed^; see below. 

Nevertheless, the basic observation of the (almost) uni- 
versal behavior inside the SG region, which we find when 
looking at the GS energy, becomes even more apparent, 
when studying the behavior of the DW energy, which we 
do in the next section. 



B. Domain- wall energy 

For all samples considered in the previous section, 
we have also calculated the DW energy as defined in 
Eqn. dH). In Fig. [3 the average value {E^w) of the DW 
energy is shown as a function of system size for selected 
values oi p. Close to the phase transition pc = 0.103, the 
systems exhibit a high degree of ferromagnetic order for 
small sizes, leading to relatively large values of the av- 
erage DW energy. This explains why very large system 
sizes were needed in Ref. [ij to determined the location 
Pc = 0.103 of the phase transition precisely. For larger 
values p > 0.2, the ferromagnetic correlation length is 
small. 

In Fig.[Slthe average absolute value {\Edw\) of the DW 
energy is shown as a function of system size for selected 
values of p. For very small values of p, close to the phase 
transition, this DW energy increases with system size for 
small system sizes. Hence, if only slow algorithm were 
available, the system would look like exhibiting an or- 
dered phase {i.e. 6 > 0). When going to larger system 
sizes, it becomes apparent that this is only a finite-size 
effect. For intermediate values of p, the DW energy de- 
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FIG. 5: The average value of the DW energy as function of 
system size for selected values of p. Error bars are very small 
for values (Edw) > 0.1 but relatively large where the DW 
energy is small. Thus, they are omitted for readability. Note 
the double-logarithmic axes. Lines are guides to the eye. 



creases with growing system size, but no saturation of 
{\Edw\) is visible on the accessible length scales, hence 
the data look similar to the results for Gaussian systems^ 
with 6 < 0. When looking at the results for p > 0.25, 
it becomes apparent that converges to a finite 

value. This is due to the quantized nature^ of the pos- 
sible values for Edw Therefore, it is clear that also for 
intermediate values of p, the DW energy converges to 
a finite plateau-value as well, but for much larger sys- 
tem sizes. Hence, the model provides a striking example 
that finite-size corrections can persist up to huge length 
scales. Note that the convergence to a non-zero plateau 
value is compatible with the results from the previous 
section, where we have found that inside the SG phase 
the scaling function Eqn. ([7]) with 9 — describes the 
data well everywhere. 

Due to the quantization of the possible DW energies, 
a value (|-Edu;|) close to zero means that many system 
exhibit actual zero DW energy Edw = 0. We have 
also studied the finite-size dependence of the probabil- 
ity P{Edw = 0) for a zero-energy DW as a function of 
the system size, as shown in Fig. [71 One observes that 
zero-energy domain walls are very common. Their prob- 
ability of occurrence increases with the concentration p 
of the antiferromagnetic bonds (until p — 0.5 due to the 
symmetry of the model in p, (1 — p)) and with the sys- 
tem size leading to a convergence to a limiting value for 
L CO. 

The simplest assumption for the convergence of {\Edw\) 
and P{Edw = 0), at given values of p, to their limiting 
values for i — > oo is a power law 
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FIG. 6: The average absolute value of the DW energy as 
a function of system size for selected values of p. Note the 
double-logarithmic axes. Lines are guides to the eye. 
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FIG. 7: The probability P{Edw = 0) for a zero-energy domain 
wall as function of system size for selected values of p. All 
error bars are at most of the symbol size. Lines are guides to 
the eye. 



The exponent uj* describes the rate of convergence. It 
depends on the fraction p of the antiferromagnetic bonds 
and on the aspect ratio R, as well as the constants A* 
and B* . We have included the dependence on R because 
below we also study values of i? < 1 , but for the moment 
we stick to i? = 1 . 

When fitting Eqn. ^ to the DW energy data for 
p > 0.15, i.e. far enough away from pc, and for sizes 
L > 10, we obtain results as shown in Fig. [51 For small 
concentrations p, the exponent is less negative than 
for p close to 0.5, corresponding to needing larger sizes 
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FIG. 8: The results of fitting Eqn. (|9} to the data {R = 1): 
exponent ui* (main plot) and limiting value A* (inset). 



to reach the asymptotic value A* . On the other hand A* 
itself seems to depend only weakly on p, confirming the 
notion that everywhere inside the SG region the DW en- 
ergy reaches a non-zero plateau value. The dependence 
of Lu* appears like a strong violation of universality inside 
the SG region. To investigate this assumption, whether 
it is a true non-universality or just a finite-size scaling 
effect, we have performed fits which also include correc- 
tions to scaling. We do this by fitting the data to a form 



F{R,p,L) = A{R,p) + BrniR,p)L- 



(10) 



m— 1 



F{R,p,L) represents any function of P{Edw)- Two ex- 
amples that we chose to study were (|i?c(u)|) and the prob- 
ability P{Edn,^0). 

The number of correction-to-scaling terms included, 
m*, was chosen according to the amount of data to be fit. 
Here we have tried 2 and 3. Of course, the values of the 
coefficients A and Bm which are found by the fits depend 
somewhat on the choice of m* . The computed statistical 
errors do not include any allowance for systematic errors 
due to the choice of the form of the scaling function. For- 
tunately, it turned out that the value of A is relatively 
insensitive to the choice of m* . 

In general, whether we chose m* to be 2 or 3, the 
computed Bm coefficients did not all have the same sign. 
This effect is caused by the behavior at small L. For this 
reason, there is no well defined prescription for deciding 
what the best choice of the exponent lo is. Wo have not 
included w in the set of free fitting parameters. Since it 
is the purpose of this section to show that the behavior 
inside the full SG region is universal, we have selected a 
value for w, such that the data for all values of p can be 
fit. Thus the strong finite-size effects close to Pc can be 



explained also. We have chosen the value lu = 2/3, which 
is similar to the values used previously in other work 
on this modelf^i^ empirically. For a critical point with 
Tc > 0, the general theory of finite-size scaling'^* requires 
that LU = 1/v, where v is the correlation length exponent. 
That result does not apply here, however, since Tc = 0. 
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FIG. 9: (color online) Finite-size scaling fits for p — 0.5 and 
R=l: (a) {\Ed^\) vs. L"^/^; (b) probability of Ea^ = vs. 
rpj^g error bars show one standard deviation. 

As an example, our finite-size scaling fits (using m* = 
3) to the data using Eqn. (fTO]) for the case p — 0.5 and 
R — 1 are shown in Fig. [SI the data for all system sizes 
was included in the fits. In the figures, we show F{R, L) 
as a function of L~^/^. One observes a straight line for 
L~^/'^ 0, showing that the choice of the exponent is 
compatible with the data. Note that the same R = 1 
data was originally^ fit (only for the case p = 0.5) by 
making different assumptions which are not consistent 
with Eqn. (fTO]) . 

Fig. [TUf a) shows the fit coefficients for {\Edw\): as de- 
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FIG. 10: (color online) Fit coefficients (defined by Eqn. [TO)) 
for 7? = 1: (a) {\Ea^\) vs. p; (b) P{Ea^ = 0) vs. p. 



fined by Eqn. pH)) . as a function of p for R = 1. The 
quality of fit was good everywhere, proving that the ap- 
parent non-universality visible in the dependence of the 
effective exponent uj* is due to the presence of corrections 
to scaling. We again see that the value of A, the estimate 
for L ^ oo, changes very little over this range of p; the 
values are for p > 0.3 very similar to the fit according to 
Eqn. ©J but they are larger in the present fit for small 
concentrations, p < 0.3. We also note that there are sub- 
stantial variations in the Bm coefficients. In particular, 
alternates in sign for the smaller values of p. In Fig. 
[TUfb) the results of fitting P{Edw — 0), which shows the 
same type of behavior, are given. In this case the signs 
of i?2 and B^ alternate. 



FIG. 11: Domain walls of subsystems are not independent: 
Comparison of a minimum-energy DW for a system of size 
100 X 200 (left) having Edw = 2 and of the same system cut 
into two systems of size 100 x 100 (right), having Edw = 0. 
Note that in the left figure, the DW wraps around the system 
in the a;-direction. Note also that the DWs are highly degen- 
erate and the algorithm does not find DWs in a controlled 
way. A "jump" in the DW at the right is visible. This figure 
illustrates what a typical DW looks like. The important point 
here is the difference in the DW energies. 



C. Aspect-ratio scaling 

The main assumption when using AR scaling Eqn. ^ 
in the limit R —^ is that different subsystems of size 
M X M are independent of each other, which leads to 
a Gaussian distribution of the DW energies Edw in the 
limit R ^ 0. Nevertheless, one can easily imagine that 
the DWs in different parts of the system are not truly 
independent of each other. This is demonstrated for a 
sample system in Fig. [TTJ where the DW energy of each 



subsystem (l)/(2) of size 100 x 100 is 



(l)/(2) _ 



= 0, but 



the DW energy of the full system is Edw = 2. Also, in a 
previous studj*^^, the validity of AR scaling could not be 
established for the ±J model. Here, we will show that 
AR scaling indeed worlts also for this case, i.e. the dis- 
tribution of the DW energies becomes indeed Gaussian, 
despite the non-independencies of the domain walls in- 
side the different subsystems. Because the shape of the 
P{Edw) distribution is changing with R, we expect that 
the AR scaling will not be perfect in the range of our 
data. As we shall see, however, the deviations from the 
Gaussian fit become smaller than our statistical errors 
for R < 1/4. Thus we are able to verify that we are 
approaching the predicted scaling limit. 

First, we concentrate onp = 1/2 where the distribution 
P{Jij) of Eqn. ([3]) is symmetric about zero. Thus in 
this case the distribution P{Edw) is (ignoring statistical 
fluctuations) also symmetric around zero, for any values 
of L and M. Therefore, we used p = 1/2 to collect data 
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for sets of random lattices with sequences of L and M 
having the aspect ratios R = 1, 1/2, 1/4, 1/8, 1/16 and 
1/32. For each value of R, we then used finite-size scaling 
to extrapolate to the large lattice limit. 




FIG. 12: The symbols display the distribution of DW energies 
for the case L = M = 40 (i. e. _R = 1) and p — 0.5, the lines 
shows the result of a fit to a Gaussian. Strong deviations from 
a Gaussian distribution are visible at P{Edw ~ 0) (main plot) 
and in the "tail" of the distribution (inset). 

We start by studying L x L samples {R ~ 1), to es- 
tablish that in this case the distribution of the DW en- 
ergies is indeed not Gaussian. (Otherwise, it would not 
be surprising if it is also Gaussian for R ^ 1.) The 
distribution for L = 40 is shown in Fig. [T^ This dis- 
tribution is obtained from an average over 39000 real- 
izations. We have fit^ a Gaussian to the data. (Note 
that the fitting procedure ignores the fact that the en- 
ergies are quantized.) Strong deviations from Gaussian 
behavior are visible at P{Edw = 0) (main plot) and in 
the "tail" of the distribution (inset). Correspondingly 
the quality of the fit, i.e. per degree of freedom, 
as given by our fitting program^ in this case, is very 
high: x^/ndf= 78. Even larger deviations are expected 
in the tail. To observe these, more sophisticated tech- 
niques would be needed^i^, which is beyond the scope 
of this work. The strong deviation from Gaussian behav- 
ior is not a finite-size effect, as may be seen from Fig. 1131 
where the circles display x^/ndf as a function of L for 
the R = l case. 

The reader should also note that we find no evidence 
for any "sawtooth" structure in the data. This is in 
marked contrast to the case of periodic or antiperiodic 
boundary conditions in the y direction, which causes al- 
ternating values of Edw to have zero probability. There- 
fore P{Edw) does not become independent of the bound- 
ary conditions even in the limit of large lattices. This 
is an aspect of the topological long-range order— which 
exists in this model when the boundary conditions are 
periodic along both x and y. 
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FIG. 13: Quality of the Gaussian fit for the DW energy dis- 
tribution measured by x^/ndf as a function of the block size 
L for p = 0.5 (circles: R = 1, squares: R = 1/32). 
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FIG. 14: The symbols show the distribution of DW energies 
for the case L = 40, M = 1280 {i.e. R = 1/32) and p = 0.5, 
the lines shows the result of fitting to a Gaussian. Main plot: 
linear ordinate, inset: logarithmic ordinate scale. The data 
match the Gaussian very well. 

In Fig.[T3]we show the probability distribution P{Edw) 
for a set of lattices with p = 1/2, R = 1/32 and L = 40. 
In this case, there are 90,000 random samples in the data 
set. Again, the distribution was fit to a Gaussian^. This 
fit is quite good, as demonstrated by the fact that the 
value of x^/ndf are is close to one. Note that the precise 
value of x^/ndf depends on the number of points in the 
tail of the distribution that are included in the fit, which 
is somewhat arbitrary. The Gaussian behavior here is 
not a finite-size effect either. This can be seen from Fig. 
I13[ where the square symbols display x^/ndf as function 
of L for the R— 1/32 case. 

When p = 0.5 the configuration averages of all the odd 
moments of F(Edw) must vanish. Then the dominant 
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FIG. 15: (color online) Kurtosis vs. L 
is scaled logarithmically. 
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TABLE II: Probability of Edw = for p = 0.5, extrapolated 
to L = CO, as a function of the aspect ratio, R. The m* and 
A are defined by Eqn. (|10|) . and the eA axe statistical error 



estimates. 


R 


m* 


A 


eA 


A{R)/A{2R) 


1.0 


3 


0.55921 


0.00294 




0.5 


3 


0.33215 


0.00261 


0.594 


0.25 


2 


0.21619 


0.00135 


0.651 


0.125 


2 


0.14736 


0.00139 


0.682 


0.0625 


2 


0.10264 


0.00144 


0.697 


0.03125 


1 


0.07053 


0.00069 


0.687 



Again, for R < 1/4 the agreement is about as good as 
could be expected. 

We have also studied the AR approach for p — 0.15 
and found similar results (not shown). For aspect ratios 
R < 1/4 the assumptions of a Gaussian distribution of 
the DW energies is again valid. Hence, we believe that 
inside the entire SG region the behavior is consistent with 
the assumptions of the AR approach. 



TABLE 1: (l-Bdml) for p — 0.5, extrapolated to L = oo, as a 
function of the aspect ratio, R. The m* and A are defined by 
Eqn. (|10|l . and the eA are statistical error estimates. 



R 


m* 


A 


eA 


A{R)/A{2R) 


1.0 


3 


0.96025 


0.00792 




0.5 


2 


1.81967 


0.00471 


1.895 


0.25 


2 


2.89085 


0.00671 


1.589 


0.125 


2 


4.28859 


0.01013 


1.484 


0.0625 


2 


6.21247 


0.02387 


1.449 


0.03125 


1 


8.98954 


0.03402 


1.447 



contribution to non-Gaussian behavior will be the kur- 
tosis, which simplifies to 7714/(7/12)^ — 3, where rrii is the 
i*'' moment of P{Edw)i under these conditions. Fig. fTSl 
shows the behavior of the kurtosis of the P{Edw) distri- 
butions at p = 0.5 as a function of L and R. It shows 
that for i? < 1/4 the kurtosis is negligible, except when 
L < 8, quantifying the convergence towards a Gaussian 
distribution for small R. 

This convergence can be also seen when fitting {\Edw\) 
to Ec^n. pTI)l . now for aspect ratios R < \. In Table I 
we show the results of the finite-size scaling fits at p = 
0.5. According to Eqn. ([S]), the prediction of AR scaling 
theory is that the ratio A{R)/A{2R) should approach 
\/2 « 1.414 for small R. Considering the statistical errors 
and the corrections to scaling, the agreement is quite 
satisfactory. Here we observe that for i? < 1/4 the scaling 
assumption is reasonably accurate. 

Table II shows the results of the finite-size scaling fits 
for the probability that Edw — Q sX p = 0.5. Since, from 
the central limit theorem, for the sum Y, oi K suitably iid 
integers, P(S = 0) if via K ^ R^^ the prediction 
of AR scaling theory for this quantity is that the ratio 
A{R)/A{2R) should approach 1/^/2 w 0.707 for small R. 



D. Discussion of AR scaling 

At the end of Section II we argued that for small R the 
system may be thought of as a set of L x i subsystems 
stacked in the y direction. This implies that the behav- 
ior of Edw in this limit should obey Eqn. ([6|) with — Q. 
As we have seen, our numerical results are indeed con- 
sistent with this argument. However, the hand-waving 
argument does not specify how the subsystems are to 
be connected to each other. Before actually seeing the 
numerical results, we were far from confident that this 
argument would correctly predict what was found. 

There is another apparently reasonable argument, 
which leads one to expect that the results for small R 
should not give this answer. As i? ^ 0, we know that the 
width of the Edw distribution will diverge. This means 
that for most of our sample lattices \Edw\ ^ 1 in this 
limit. One might have thought, therefore, that the fact 
that Edw is quantized would no longer matter. Thus one 
might have guessed that the system would behave in this 
limit as if 6* « —0.28, the value for an unquantized dis- 
tribution. 

Therefore, our results indicate that a Bohr correspon- 
dence principle, that in the limit of large quantum num- 
bers the results should be indistinguishable from those of 
an unquantized system, does not apply. How can we un- 
derstand this? In our opinion, what we learn from this 
is that = is directly connected to the existence of 
topological long-range order— at T = 0. 

One might try to argue that R = 1/32 is not small 
enough, and that for even smaller R one would indeed 
find a crossover io k, —0.28. Such a viewpoint is sug- 
gested by the recent work of Jorg et ali^ On the other 
hand a more detailed analysis^^ of the low temperature 
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behavior of the specific heat of the ±J model does not 
support this interpretation. 

IV. SUMMARY 

We have studied the statistics and finite-size scahng be- 
havior of GS and DW energies for the 2D Ising random- 
bond spin glass with an mixture of +1 and —1 bonds, 
where p is the concentration of antiferromagnetic bonds. 
By using sophisticated matching algorithms, we can cal- 
culate the GS energy and DW energies for quite large 
systems exactly, which allows us to study the scaling be- 
havior very precisely. We find that the scaling behavior 
of the GS energy can be described everywhere inside the 
SG region by the same simple scaling function Eqn. ([7]), 
as predicted by Ref. [H. Furthermore, when looking at 
the fluctuations of the GS energy, we find, after care- 
fully taking into account finite-size corrections, that they 
follow the predicted simple i~^-scaling with an ampli- 
tude which has a cusp singularity at the ferromagnetic- 
SG transition pc- The singularity is described by an ex- 
ponent which is close to 2/3. 

The behavior of the DW energy is characterized close 
to Pc, by huge finite-size effects. Nevertheless, the results 
are compatible with a convergence of (jiSdiol) to a finite 



plateau value everywhere inside the SG phase. This con- 
vergence can be described, again universally inside the 
SG phase, by a single exponent w ~ 2/3, just by taking 
into account higher-order corrections to scaling. Finally, 
we have studied L x M rectangular lattices with aspect 
ratio R between 1/32 and 1. We have demonstrated by 
an example that one assumption underlying AR scaling, 
i.e. the assumed independence of the DW energies of dif- 
ferent blocks, is not strictly valid. Nevertheless, we find 
that for large lattices the probability distribution of E^iyj 
in the SG region of the phase diagram approaches for 
R < 1/4 a Gaussian centered at Edw = 0. Hence, in the 
small R limit the behavior obeys the AR scaling predic- 
tions of Carter, Bray and Moorej^ 
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